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Abstract 


We have evaluated the centrifugal force acting on a fluid element and the ellipticity 
of the fluid configuration, which is slowly rotating, using the Hartle-Thorne solution 
for different equations of state. The centrifugal force shows a maximum in every case, 
whereas the reversal in sign could be seen in only one case, and the system becomes 
unstable in other cases. The ellipticity as calculated from the usual definition shows 
maxima, whereas the dehnition obtained from the equilibration of the inertial forces, 
shows a negative behaviour, indicating that the system is prolate and not oblate. 
This prolate shape of the configuration is similar to the one earlier found by Pfister 
and Braun for a rotating shell of matter, using the correct centrifugal force expression 
for the interior. The location of the centrifugal maxima gets farther away from the 
Schwarzschild radius (R^) as the equation of state gets softer. 


Introduction 

One of the most important aspects in the discussion of the dynamics of a rotating 
neutron star is the equation of state of matter contained within. There have been 



a large number of studies in this regard starting from the pioneering paper of Op- 
penheimer and Volkoff |]I| in 1939 and later by Pandharipande (1971) p| and several 
others over the last twenty hve years. In spite of all the efforts put in, it is still enig¬ 
matic as to the correct nature of the equation of state and its transitions within, as it 
involves the actual short range interactions about which there exists no exact theory 
as yet. However, there have been several different equations of state, considering 
the nuclear interactions at various levels and potentials, and it is important to see 
the effects that they may have in the actual structure and dynamics of a neutron star. 


It is indeed well known that a rotating body breaks spherical symmetry and the 
equilibrium configuration would be a Maclaurin spheroid or a Jacobi ellipsoid due to 
ellipticity arising from the anisotropy of the equilibrating forces. Our interest is to 
consider the dynamics of the interior of ultra compact bodies in the realm of general 
relativity as the gravitational potential is quite large. Though the original formula¬ 
tion of general relativity does not express the equilibrium in terms of forces, a recent 
3-1-1 formulation allows one to introduce the language of forces without affecting any 
of the relativistic effects. On the other hand, through the language of forces one 
gets some new results which were not apparent earlier. One of the interesting new 
features is the realisation that the centrifugal force acting on a fluid element of the 
conhguration has a maximum value and that it also reverses its direction at a value 
of r, outside the event horizon |P], [Q. Further, it has also been seen that the el- 
lipticity of the spheroidal conhguration attains a maximum [Q, a result which was 
earlier shown by Miller et ah |p. As the change in ellipticity would have to do with 
the shape of the compact object, it is indeed interesting to check how the ellipticity 
and centrifugal force change with different equations of state. As shown earlier in 
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the case of homogeneous distribution, the centrifugal force reverses for a conhgura- 
tion with R/Rg ~ 1.45, and the ellipticity for the same conhguration is maximum at 
R/Rs ^ 2.75. 


In this, we extend this study to the case of an inhomogeneous density distribution 
with four different equations of state as given below: 


(A) Pandharipande (hyperonic matter) p, 

(B) Wiringa, Fiks and Fabrocini’s beta-stable model: UV14-I-UVII |]^, 

(C) Walecka’s model corresponding to pure neutron matter [§, and 

(D) Sahu, Basu and Datta’s model based on the chiral sigma model p. 

Model (A) is Pandharipande’s (hyperonic matter) EOS which studies the behaviour 
of dense matter using a many body theory based upon the variational approach. Hy- 
perons are considered one of the baryonic constituents of the neutron star’s interior. 
Model (B) considers the main constituents of neutron star matter to be neutrons, pro¬ 
tons, electrons and muons. This is one of the three models given by Wiringa, Fiks and 
Fabrocini which includes three-nucleon interactions using a non-relativistic approach 
based on the variational method. Walecka’s model (C), which characterizes the effec¬ 
tive interaction by the meson parameters (masses and the coupling constants), is for 
pure neutron matter and emphasises the scalar and vector meson exchange interac¬ 
tions by considering the relativistic approach. Model (D) is a held theoretical EOS 
for neutron rich matter in beta equilibrium based on the chiral sigma model. 
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Model A is the softest equation of state among the considered EOS, whereas the most 
stiff EOS is model D. B and C are intermediate EOS of which C is the stiffer of the two. 


Using the ACL formalism of 3+1 conformal slicing of the space time 

= df + gooidt + 2ujadx°‘Y, ( 1 ) 

with dP representing the positive definite metric of the absolute 3-space g^^dx^dx" 
(where Greek indices take values from 1 to 3 and'denotes the quantities in absolute 
3-space), the four force acting on a fluid element of a slowly rotating perfect fluid 
distribution, represented by the spacetime metric 

ds^ = -e^'^dt^ + , ( 2 ) 


in a locally non rotating frame (LNRF), is given by 0 


/o = <h + p)U>^df,Uo + 


( 3 ) 


/„ = $ (p + p) 


' ~ ~ ~ 

2<h 


+ Kp,, 


( 4 ) 


Here <h represents the gravitational potential {—goo), P the pressure, p the energy 
density and the velocity four vector of a fluid element. An overhead tilde rep¬ 
resents the quantity in the projected 3-space of the optical reference geometry. The 
terms a and {Mq/ 2^)da^ represent the centrifugal (E),/) and gravitational 

acceleration (Eg/) respectively and have the form 


Fcf 

F, 






( 5 ) 


for metric (2). 
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The Newtonian force balance equation 


dequator 9polei,^ C 

may be expressed in general as 


- F,f = F,,(l - 


( 6 ) 

(7) 


with 'Fge and 'Fg^ representing the accelerations at equator and pole respectively. 
In the limit of slow rotation, the ellipticity is given by 

1 o 


e = -e 

2 


( 8 ) 


where e is the eccentricity of the spheroid. Using equations (7) and (8), ellipticity is 
seen to be 


For the Hartle-Thorne 


1 

e = - 
2 


10 metric 




Fcf = rV(l/r - z/'/2), 

Fge = ~ ^ 2 ) + 2hQ — 

^gp ~ -e^°[z/Q(l + 2 ho + 2/12) + 2 hQ + 2/2,2], 

yielding for the ellipticity function: 

2 — 2 

e = 3(ft2 + {2/rv'„ - 1). 


(9) 


( 10 ) 

( 11 ) 

( 12 ) 


(13) 
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Results and Discussions : 


Table 1 shows the location of extrema for the centrifugal force as well as for the el- 
lipticity and their respective values for both the Hartle-Thorne definition {Ih-t) and 
our dehnition (e). Figures (l)-(3) give the plots of centrifugal force Fcf and the two 
ellipticities (expressed in dimensionless units and respectively) as a 

function of R/Rg, Rs being the Schwarzschild radius. 


As is seen, the centrifugal force keeps increasing as the conhguration size gets smaller 
and then attains a maximum somewhere between R/Rg = 2.1 and 2.3, for different 
equations of state. However, unlike in the case of a homogeneous spheroid where the 
centrifugal force reverses sign at i? = lA5Rg, with these different equations of state 
the reversal is seen only in the case of Wiringa et ah model (B), at i? ~ 1.454 Rg. 
In other cases the equilibrium configuration becomes unstable before reaching the 
value R = 1.5Rg, which in fact is the radius of the orbit of the particle for which the 
centrifugal force is zero in the Schwarzschild space time. 


Considering the nature of ellipticity it appears that the behaviour is different for the 
inhomogeneous distribution than the case of homogeneous distribution. As seen from 
Fig. 2, the ellipticity (e) keep reducing as the conhguration gets smaller, becomes zero 
(meaning that the shape becomes spherical) and further on gets to a negative value. 
However, the negative value attains a minimum and then again the ellipticity starts 
increasing. 


The fact that the ellipticity starts decreasing and further becomes negative, could be 
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due to the reason that the force balance equation used in dehning the ellipticity is 
perhaps true only for a homogeneous distribution of the fluid, whereas we have in the 
above varying density conhgurations. However, a point that needs to be checked care¬ 


fully is that when Phster and Braun used the correct centrifugal force expression 
for obtaining the solution for the interior of a mass shell, they found that a proper 
boundary £t of the exterior and the interior solutions for the shell, was possible only 
if the conhguration is prolate rather than oblate. Here, in our approach also we start 
from the equilibration of the forces within the framework of general relativity and get 
prolate conhguration for distributions with inhomogeneous density. 


It is also worth noting that for the same conhgurations, when ellipticity is dehned 
in terms of the radii of the object with constant surface density embedded in a 3- 
dimensional hat space, a la Hartle-Thorne [ im , one gets the ellipticity maximum as 
was in the case of a purely homogeneous conhguration. In this case, the location of 
the maxima changes with the equation of state, shifting inwards as the equation of 
state gets stiher. As the equation of state of any conhguration describes the pressure- 
density relation, the equilibrium of the conhguration in a sense depicts the balancing 
of the various forces like gravity, material binding and the centrifugal force. As the 
equations of state gets stiher the intra-nucleonic forces, which are ehectively repul¬ 
sive at very short range become larger, thus requiring the conhguration to get more 
compact before similar behaviour of ellipticity extrema is attained. However, a mat¬ 
ter of concern is the result concerning the diherence in behaviour of the ellipticity 
function in the two diherent treatments, for the inhomogeneous distribution while 
the conventional Hartle-Thorne way of dehning it via embedding in a 3-hat geometry 
shows the function to be positive, dehning it through the balancing of inertial forces 
a la Maclaurin and Newton, shows it to be negative. It is important to look deeper 
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into this question to ascertain whether the treatment of defining inertial forces for 
a fluid configuration has to be different from the approach used for a single particle 
dynamics, particularly for inhomogeneous distributions. 
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Figure Captions 


Fig. 1 Centrifugal force Fcf in units of for decreasing values of radius R in 

terms of Schwarzschild radius Rg for various equations of state. EOS A is the 
softest and EOS D represents the stiffest among the four considered equations 
of state. 

Fig. 2 Ellipticity e derived in optical reference geometry, in units of 
Fig. 3 Ellipticity en-r as dehned by Hartle and Thorne, in units of J'^/M*. 
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Table Captions 


Table 1 Location of extrema for the centrifugal force {Fcf) and ellipticity (e) (column 
7, 5) and their values (column 6, 4) for the equations of state (Models A, B, 
C, D) considered in this paper as well as for the homogeneous distribution. 
Column 8 gives the location of reversal in e. The radius is expressed in terms of 
Schwarzschild radius Rs{= 2M), whereas the ellipticity and the centrifugal force 
are expressed in the dimensionless units of and respectively. 
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Table 1 


EOS 

^H-T 

TD 

■^en-T 

e 

R-e 

Fcf 

Rp f 

Re 

^rev. 

A 

0.953 

3.278 

-0.257 

1.603 

0.0335 

2.234 

2.297 

B 

0.849 

2.948 

-0.403 

1.610 

0.0271 

2.204 

3.301 

C 

0.832 

2.857 

-0.395 

1.769 

0.0261 

2.180 

3.631 

D 

0.813 

2.625 

-0.388 

2.014 

0.0252 

2.112 

4.343 

Homo. 

0.761 

2.3 

1.207 

2.75 

0.0157 

2.1 
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